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The Nambu-Jona-Lasinio jet model is extended to accommodate hadronization of a transversely 
polarized quark in order to explore the Collins effect within a multihadron emission framework. 
This is accomplished by calculating the polarized quark spin flip probabilities after a pseudoscalar 
hadron emission and the elementary Collins functions. The model is then used to calculate the 
number densities of the hadrons produced in the polarized quark's decay chain. The full Collins 
fragmentation function is extracted from the sine modulation of the polarized number densities 
with respect to the polar angle between the initial quark's spin and hadron's transverse momentum. 
Two cases are studied here. First, a toy model for elementary Collins function is used to study the 
features of the transversely polarized quark-jet model. Second, a full model calculation of transverse 
momentum dependent pion and kaon Collins functions is presented. The remarkable feature of our 
model is that the 1/2 moments of the favored Collins fragmentation functions are positive and peak 
at large values of z but decrease and oscillate at small values of z. The 1/2 moments of the unfavored 
Collins functions have comparable magnitude and opposite sign to the favored functions, vanish at 
large z and peak at small values of z. This feature is observed for both the toy model and full 
calculation and can therefore be attributed to the quark-jet picture of hadronization. Moreover, 
the transverse momentum dependencies of the model Collins functions differ significantly from the 
Gaussian form widely used in the empirical parametrizations. Finally, a naive interpretation of 
the Schafer-Teryaev sum rule is proven not to hold in our model, where the transverse momentum 
conservation is explicitly enforced. This is attributed to the sizable average transverse momentum 
of the remnant quark that needs to be accounted for to satisfy the transverse momentum sum-rule. 

PACS numbers: 13.60.Hb, 13.60.Lc, 13.87.Fh, 12.39.Ki 
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I. INTRODUCTION 

Deep inelastic scattering (DIS) has been a powerful 
tool with which to explore the structure of hadrons. Par- 
ticularly, experiments using the semi-inclusive deep in- 
elastic scattering (SIDIS) have advanced our understand- 
ing of the underlying partonic structure of nucleon in mo- 
mentum space, allowing access to both the longitudinal 
and transverse momentum distribution of quarks and glu- 
ons. Using polarized probes and/or targets has allowed 
us to extend these explorations to the spin dependence 
of the parton distributions. In particular, the naively 
time-reversal odd Sivers distribution function has been 
measured to be nonzero [TH3]. The SIDIS cross sections, 
both polarized and unpolarized, are a convolution of par- 
ton distribution functions with elementary parton-probe 
scattering amplitudes calculable using perturbative QCD 
and parton fragmentation functions c 4]. The complete 
tree level expressions at leading order have been written 
in Ref [5]. Thus our detailed understanding of nucleon 
structure from SIDIS is hinged on our knowledge of frag- 
mentation functions. Moreover, the naive T-odd distri- 
butions are convoluted with T-odd fragmentation func- 
tions, notably the Collins fragmentation function H^. 
The first direct measurements of the Collins mechanism 
have been performed by the Belle collaboration using 



hadron pair production in e + e collisions 

There has been a vast amount of work done to calculate 
the Collins function within theoretical models. Most no- 
tably, the first calculation of a nonvanishing Collins func- 
tion for pion production within a quark spectator model, 
using the interference of one-pion-loop amplitudes with 
the tree level amplitude, was preformed in Refs. [H [9j. 
The resulting Collins fragmentation function appears to 
be too small to describe the experimental data as shown 
in Ref. [TO], where the contributions from initial quark 
rescattering were largely canceled by the quark-pion ver- 
tex loop corrections. Later, a similar mechanism was 
used to calculate the Collins function, through inter- 
ference of the gluon rescattering and gauge-link correc- 
tion amplitudes interference with the tree level amplitude 
Refs. [TUrfT^] . The resulting Collins functions are in gen- 
erally good agreement with the data, but this approach 
lacks the ability to produce the unfavored fragmentation 
functions, including Collins functions. The experimental 
results from HERMES, COMPASS and JLab are strongly 
suggesting that the unfavored Collins functions have a 
similar size and an opposite sign to that of the favored 

ones mo nam]. 

The NJL-jet model of Refs. [151118] has been used 
to calculate the quark fragmentation function within 
the quark-jet hadronization picture of Field and Feyn- 
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man |19j . using the elementary quark-hadron splitting 
functions calculated within the effective quark model of 
Nambu and Jona-Lasinio (NJL) [20]. The advantage of 
the model is that it has no free parameters adjusted to 
fragmentation data and it naturally describes the unfa- 
vored fragmentation functions. The use of Monte Carlo 
(MC) methods in Refs. [T7] [TS] allows for inclusion of 
hadronic resonances and their decays as well as the trans- 
verse momentum dependence of both parton distribution 
and fragmentation functions. The model was also used to 
calculate the dihadron fragmentation functions and their 
evolution in Refs. [2T][22]. Here we propose to extend the 
NJL-jet model by including the spin of the fragmenting 
quark to calculate the Collins function, both favored and 
unfavored. 

Here we consider the fragmentation of the transversely 
polarized quark, as depicted schematically in Fig. [T] The 
goal is to calculate the T-odd Collins fragmentation func- 
tion using the NJL-jet Monte Carlo method of Ref. [18] . 
Following the "Trento Convention" [23], the fragmenta- 
tion function of the transversely polarized quark q to un- 
polarizcd hadron h carrying the light-cone momentum 
fraction z and transverse momentum P± with respect to 
quark's momentum k can be expressed as a sum of two 
terms 



D h/q r(z,Pl,^)=D'l /q (z,Pl 



(1) 



+ Hl h/q (z,Pl) y "' " * ±: 



(fe x P ± ) ■ S q 



zm h 



where k and S q are the momentum and the spin vector 
of the fragmenting quark, rn/j is the mass of the pro- 
duced hadron. The unpolarized fragmentation function 
is denoted as D ] ^ q {z,P\) and H^'^ q (z, P±) is the so- 
called Collins function. The spin dependent geometric 




FIG. 1. Illustration of the three dimensional kinematics of 
transversely polarized quark fragmentation. The fragmenting 
quark's momentum k defines the z-axis with its transverse 
polarization spin vector S q along x axis. The emitted hadron 
has momentum P with the transverse component P± with 
respect to the z-axis. The polar angle of hadron's momentum 
P with respect to the zx plane is denoted by p. 



factor multiplying the Collins function in Eq. ([I]) can be 
expressed in terms of the polar angle ip of the vector P 
using the 3- vector components, yielding for the polarized 
fragmentation function 



D h/qf (z,Pl,v) = D , { /q (z,Pl 



(2) 



H 



±h/q 



P±Sq 



zm h 



sin(y). 



The measure for the above number density is given as 



dN h/qT (z,P x ) = D h/qt (z,Pl,<p) dz d 2 P ± 



D h ^(z,Pl) 



H 



zm h 



(3) 

dP 2 

dz—^dp. 



Integrating over the azimuthal angle p we recover the 
expression for the number density of the unpolarized 
quark 



dN h/qf (z,Pl) 



dp N h/qT (z,P 2 ,ip) 



(4) 



7rD^ /q (z,Pi) dz dPi = dN h/q (z,P'l). 



On the other hand, if we perform an integration over 
Pj_, we arrive at an expression that allows us to extract 
the 1/2 moment of the Collins function H^S^i q ?(z) 



dN h/q t(z,ip) = D h / q f(z, ip) dz dip 
dPl N h/q ,(z,Pl,<p) 



(5) 



1 

'2tt 



D h /\ Z ) 



dz dip, 



where 



D\' q {z) 



n dPlD^iz^l), 
Jo 



H 



±(1/2) 
l(h/g) 



coo p 
T / dP 2 

la 



_L 



_L 



2zm,h 



H 



±h/q 



(z,p1). 



(6) 
(7) 



At last, integrating over both Pj_ and p> in Eq. ^ we 
arrive at the integrated number density 



dN h/q t (z) = 



dP\ 



2k 



dp N h/qJ (z,P 2 ,p>) (8) 
= D^ /q {z) dz. 



Our strategy here is to use the NJL-jet model to calcu- 
late the number densities N h / q t(z, P 2 , ip) and use the de- 
pendence on the azimuthal angle ip in Eq. (|3| to extract 
both D^ /q (z,P 2 ) and Hj- h/q (z,Pl), where the former 
can be compared with our earlier results in Ref. [IS] as 
a cross-check. To achieve this goal, we first expand the 
NJL-jet model and the associated Monte Carlo frame- 
work to accommodate the transverse spin of the frag- 



menting quark in Sec. [nj In Sec. Ill we calculate the 
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elementary unpolarized and Collins fragmentation func- 
tions, needed as an input to the NJL-jet model. In 
Sec. |IV[ we present the results for a simple toy model 
used as elementary Collins function to demonstrate the 
distinctive features of the model. In Sec. |Vj we discuss 
the full model results and present the conclusions and 
some final remarks in Sec. IVII 



II. 



quarks 



NJL-JET WITH TRANSVERSELY 
POLARIZED QUARK 



(9) 



where x — ^l are the eigenvalues of the spin projec- 
tion onto the x axis and t/(±z) are the Lepage-Brodsky 
spinors in hclicity basis 26J. Note that we have a plus 
sign in front of the second term with Xt as there is a 
sign error in the corresponding expressions of [24] . The 
spinors U x are both solutions of the Dirac equation and 
the eigenstates of the Lorentz-covariant spin operator: 
the Pauli-Lubanski vector 



The NJL-jet model employs the quark-jet mechanism 
to describe the hadronization process. Here we use the 
model to describe the fragmentation of the transversely 
polarized quarks to unpolarized hadrons, schematically 
depicted in Fig. [2] 

An important intricacy arises from the need to keep 
track of the transverse spin of the quark in the jet as 
it emits hadrons. We first examine the elementary pro- 
cess where a transversely polarized quark emits a sin- 
gle hadron, and we calculate the probability of the final 
quark's spin being parallel or antiparallel to the origi- 
nal quark's spin direction. In this work we only consider 
the emissions of pseudoscalar mesons. Further, we make 
a first order approximation by including in the calcula- 
tions of the quark spin flip probabilities the elementary 
hadron emissions only via the tree level diagram depicted 
in Fig. [4j thus neglecting any T-odd effects. The T-odd 
effects are considered here to be small relative to the 
included unpolarized term, though essential for generat- 
ing the elementary Collins function. This approximation 
can be easily improved on in the future by also including 
in the quark spin flip calculation the relevant diagrams 
that generate the elementary Collins function, such as 
that discussed in Sec. IIIIB1 

We use the kinematics depicted in Fig. [I] with the 
fragmenting quark's momentum defining the z axis. We 
denote the remnant quark as Q with momentum and spin 
vectors I and Sq . We use the Dirac spinors of [MJ [S3] 
to describe the wave functions of transversely polarized 



W = — ^Pk" 

sr>> = \[>f,>f], 



(10) 

(ii) 



where e^ up(T is the Levi-Civita tensor with convention 



e oi23 = +1- It is easy to check explicitly, that 



~m)U x = 0. 



The normalization of the spinors is 



U x (k, m)U x > (fc, m) 



8 XtX i2m. 



(12) 
(13) 

(14) 



Then the fragmenting quark's spinor is — U\(k, Mi), 
while the remnant quark's spinor can be described as 
a superposition of states with spin polarization parallel 
and antiparallel to the x axis: ^ ut — o-i U\(l,M2) + 
a_x U—\(l,M%), Then the relative probabilities of 
the quark spin flip and nonflip are determined by 
, I , cmt7 5, I'm ■ The spinor matrix elements are given by 



|[7 X ,(;,M 2 )7 5 ^ X (A:,M 1 )| 2 = 



(15) 



7 x,x' 



1 



s x,-x' 



ll + (M 2 - (1 - z)M 1 f 
1 - z 



and hence the spin nonflip and spin flip probabilities are 
proportional to 




FIG. 2. NJL-jet model including transverse momentum and 
quark polarization transfer. Here the orange double-lined ar- 
rows schematically indicate the spin direction of the quark in 
the decay chain. 



K| 2 ~ ll \a^\ 2 ~f y + (M a - (1 - z)M 1 f. (16) 

We note, however, that the spinors U x (k,m) become 
the eigenstates of W\ only if the momentum k has no 
transverse components: k x = k y = 0. Though the fi- 
nal state quark necessarily has a transverse momentum 
l± = (lx,ly) — ~ required by the momentum con- 
servation in the transverse plane, it is assumed to be 
small compared to the light-cone momentum. Thus, in 
each hadron emission step, we assume the fragmenting 
quark being in an eigenstate U x , and sample the spin 
flip probabilities according to (16). It will be discussed 
in Sees. IIVI and |V] that the results for the Collins func- 
tion simulations converge rapidly with the number of the 
hadron emissions in the region of z > 0.02, similar to the 
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FIG. 3. Quark elementary fragmentation kinematics, for an 
arbitrary hadron emission in the cascade chain. The z axis 
is defined by the direction of the 3-momentum of the original 
parent quark. 



case of the unpolarized fragmentation function that was 
examined in Ref. |17j . Thus the failure of this approxi- 
mation after many hadron emissions, when the light-cone 
momentum of the quark becomes small and comparable 
to the transverse component, affects the results only at 
extremely small values of z and can be neglected. 



Using the probabilities of the quark spin flip of Eq. ( 16 ) 



and the elementary polarized number densities, to be cal- 
culated in the next section, we perform a Monte Carlo 
simulation within the NJL-jet framework to calculate 
the average multiplicities of the produced pseudoscalar 
mesons (pions and kaons), N h / q t, as functions of light- 



cone momentum fraction z, transverse momentum square 
Pj_, and polar angle tp, in the coordinate system defined 
by the initial quark momentum and spin polarization. 
Schematically, the elementary hadron emission process is 
depicted in Fig. [3j where the vectors fc and k' denote the 
3-momentum of an arbitrary quark in the cascade chain 
before and after hadron emission with transverse com- 
ponents k± and k' ± , respectively. The emitted hadron's 
momentum is labeled by ph , where its transverse compo- 
nent with respect to fc and the z axis is denoted by p± 
and _Pj_, respectively. P± is obtained using the relation 
P± = p± + z fcj_ . The recoil transverse momentum of 
the final quark, k' ± , is calculated from momentum con- 
servation in the transverse plane, namely, 



fc 



fc', 



(17) 



In each hadron emission step we randomly sample the 
type, the relative light-cone momentum fraction z and 
transverse momentum and polar angle ip (relative to the 
spin vector of the fragmenting quark) for the produced 
hadron according to the relevant elementary number den- 
sity. We recover and record these quantities in the initial 
quark's frame, as well as sample the spin direction of the 
final quark according to (16). Then we can extract the 
number densities corresponding to Eq. ^ in the usual 
way OUCH] 



AP 2 i \ 

D h/qt (z,P 2 ,p) Az — i Ap= (N> q \(z,z + Az;P 2 ,P 2 +AP 2 ;<p lV > + Av)) 



(18) 



E 



N S i 



N^(z,z + Az,Pl,Pl+APl;ip,ip + Aip) 



Sims 



where iV^ t is the number of the hadrons h produced by the quark q that have momentum components laying within 
the regions specified in its arguments and Nsims is the number of quark decay chain simulations performed. 

III. ELEMENTARY FRAGMENTATIONS 

The quark fragmentation functions are defined using the quark-quark correlators [H [TTJ [27] 



a(z,p ± ) = y^ J %0jr eik ' 6 (o\^ T [^T^T;-^-]u-[-^-,r^T}m\h,x) 



(19) 



x (h,X\ ^(0) Zr[(r,-oo-;0r] U T [0 T , co T ; -oo~]\0) | $+=0 , 



where tp is the quark wave function, and Li^~' T ^ repre- 
sents the gauge links along the minus and transverse 
directions in light-cone coordinates. The details of the 
gauge link structure can be found in Ref. [27] , 

In this section we aim to calculate the elementary un- 
polarized and Collins fragmentation functions (or split- 
ting functions), corresponding to only a single hadron 



emission truncation of the relation ( |19| ). We use the 
framework of the NJL effective chiral quark model for 
this purpose. The corresponding unpolarized splitting 
functions, have been calculated in Refs. [15Ul7j . where 
the model parameters were also fixed. In the next sub- 
section we briefly summarize their results and refer the 
reader to the original articles for the details. For the 
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FIG. 4. Feynman diagram describing the elementary quark 
to hadron fragmentation functions. 



Collins splitting functions, we adopt the gauge link cou- 
pling mechanism to the final quark of Refs. [TUl412j . but 
use the NJL model formalism and parameters. 



A. Unpolarized splitting function 

The unpolarized fragmentation function, D^ q '{z, P± 
is defined as a trace of the quark-quark correlator 



A(z,P|)=Tr[A[z,Pj 7 +]/2. 



(20) 



The single hadron emission cut diagram for a quark 
q, to emit a meson h, carrying light-cone momentum 
fraction z, and transverse momentum p±, is depicted in 
Fig. [4] In the frame where the fragmenting quark has zero 
transverse momentum, but a nonzero transverse momen- 
tum component —p±/z with respect to the direction of 
the produced hadron [TI>1 [55] , the unregularized elemen- 
tary transverse-momentum-dependent (TMD) fragmen- 
tation functions to pseudoscalar mesons are given by 



d\ /q {z,p\) = T^ShqQ z 



16tt 3 



(21) 



V \ + [{z - l)Afi + M 2 y 



[p\ + z(z - 1)M 2 + zMl + (1 - z)m\] 

where M\ and M 2 denote the masses of fragmenting and 
remnant quarks. Quark flavor is also indicated by the 
subscripts q and Q, where a meson of type h has the 
quark flavor structure h — qQ and denotes the me- 
son mass. The corresponding isospin factor and quark- 
meson coupling constant are labeled by and ghqQ, 
respectively, and have been determined within the NJL 
model QI1CEE]. 

In the model, the transverse momentum integrated 
fragmentation of Eqs. ( 6|7 | (as well as the quark-hadron 
couplings, etc.) require an introduction of a regulariza- 
tion scheme. In this work, the integrals with two particles 
in the loop with masses /ii and fj. 2 , carrying light-cone 
momentum fractions x and 1 — x, respectively, are reg- 
ularized using the modified Lepage-Brodsky transverse 
momentum cutoff functions, as described in [17] 



[i + {Mi 2 /k\ 2 y 



(22) 



where M 2 2 is their invariant mass squared. In the frame 
where the total transverse momentum is zero, Ml 2 = 
(l4 + P±)/x + + PD/i 1 - x). The cutoff, A 12 , is 
determined by 



M? 2 < A 2 2 = 



A§ + Mi 



i4 



(23) 



where A3 is the 3-momentum cutoff, that is fixed in 
the usual way by reproducing the experimental pion de- 
cay constant. We use a light constituent quark mass 
M = 0.4 GeV, for which the corresponding 3-momentum 
cutoff is A 3 = 0.77 GeV. The strange constituent quark 
mass of M s — 0.59 GeV is obtained by reproducing 
the experimental kaon mass. The model values for the 
quark-meson coupling constants are determined from the 
residue at the pole in the quark- ant iquark t matrix to be 
g^qQ = 4.24 and g Kq Q = 4.52. 



B. Elementary Collins function 

The Collins function is defined via the trace of the 
quark-quark correlator 



£ f fc Tj 

m h 



Ht(z,Pl) = Tr[A(z,P 



'±)l(T l 75] 



(24) 



where mh is the mass of the produced hadron, = 
-ef = 1, and ei. 1 = e|? = . 

It was shown explicitly, that the simple tree level ap- 
proximation for the correlator, used to calculate the un- 
polarized quark splitting function, yields a zero result for 
the Collins function [8j[T0]. Thus, here we generate the 
model Collins function through the interference between 
the tree level hadron emission amplitude and the one with 
a gauge- link interaction with the final state quark via sin- 
gle gluon exchange, as described in Ref. [11]. Note, that 
we omit the gluon-quark rescattering and gluon loop cor- 
rection to the quark-hadron vertex, as those nearly can- 
cel each other. Also, we note that the amplitudes written 
in the spectator formalism coincide with those in NJL, 
except for the different quark-hadron coupling constant, 
the remnant " spectator" being treated as a quark in NJL 
and the loop integral regularization scheme employed in 
the NJL-jet. Thus we can simply use the results from 
Ref. [TT] with just a few modifications. The relevant di- 
agram for the NJL-jet model is depicted in Fig. [5] 

The result for the elementary Collins function is 



2a s 



s n °<? 2 m h 



1 



(27r) 4 ~ r 2 » mqLJ \ ~ z k 2 - Ml 2p\ 
x {-I 34g [(2z - 1)M X + M 2 ] 
+ I 2g \lzM x {k 2 - M\ + m\{\ - 2/z)) 
+ (M 2 - M x )((2z - l)fc 2 - m 2 h 
+M 2 2 - 2zM x (M x + M 2 ))] } , (25) 

where a s is the strong coupling constant, Cp = 4/3 is 
the cubic Casimir invariant of the color SU(3) group and 
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FIG. 5. Cut diagram describing the elementary quark to 
hadron Collins fragmentation functions. 



the initial quark's total momentum squared, k 2 , can be 
expressed as k 2 = p\/{z{l — z)) + M|/(l — z) + rri^/z. 
Since only the imaginary part of the amplitudes con- 
tribute to the above expression, the loop integrals for 
those can be calculated without the need for explicit reg- 
ularization 



2g 



2y/X(M 2 ,m h 
hig = 7t In 



In 



k 2 + M 2 - m| - y/A(M 2 , m h ) 
k 2 + Mf - m 2 + y/X(M 2 ,m h ) 



'k 2 (l-z) 

A(mi,m 2 ) = (fc 2 - (mi + m 2 ) 2 )(k 2 - (mi - m 2 ) 2 ). 



(26) 



Then the elementary polarized fragmentation function 



is given as 



,h/g 



H 



±h/g 



) P±Sq 

zm h 



sin(^), 
(27) 



where we use the multiplicative regulator of Eq. ( 22 ) for 
integrals over p 2 ^. The only remaining parameter is the 
strong coupling in Eq. (25). We consider it as a model 



parameter, that we fix to the largest value of a s — 0.444, 
that still allows for the positivity bound to be satisfied, 
namely d h j q y ^ 0. At the next-to-leading order, this 
value corresponds to a typical hadronic scale of Q 2 = 
1 GeV 2 , which is much higher than the one typically 
used as the NJL-jet model scale, namely Qq = 0.2 GeV 2 . 
Such a discrepancy between the scales for the unpolar- 
ized and Collins functions is rooted in the model for the 
Collins function employed here, which would violate the 
positivity bound if calculated at the typical model scale. 
A similar issue was encountered in the original work of 
Ref. [H], where the value of a s — 0.2 was chosen, much 
smaller than that for the scale of their model set at 
Q 2 = 0.4 GeV 2 . A completely consistent determination 
of a single scale for the polarized fragmentation function 
must involve the QCD evolution of both the TMD un- 
polarized and Collins functions. This is not possible at 
present, because the evolution equation for the Collins 
function is unknown. 

Within the NJL-jet model, the elementary split- 
ting functions are renormalized such that quark's total 
probability of emitting a hadron in each step is one: 




(a) 




(b) 



FIG. 6. Elementary renormalized unpolarized fragmentation 
function, di, (a) and Collins function 1/2 moment, 2H^ i - 1 ^ 2 \ 
(b) used in the full model calculations of Sec. [V] 



J2h f d z dp±/2 dip d h / q r(z,p 2 l , ip) = 1, the sum is over 
all hadrons the quark of given flavor can emit directly. 
These renormalized splittings will be used in the next two 
sections as input to the Monte Carlo simulations of the 
quark-jet hadronization process. The integrated renor- 
malized elementary fragmentation functions for the full 
model calculations of Sec. [V] are depicted in Fig. [6] 



IV. THE QUARK-JET EFFECTS ON COLLINS 
FUNCTION USING A TOY MODEL 

In this section we employ a toy model for the elemen- 
tary Collins function to explore the general features of the 
NJL-jet model extended to transversely polarized quark 
fragmentation, as described in Sec. JTTJ In this toy model 



we assume that 



-Ht h/q (z >P i: 



p±b q 



= 0.1 



d'; /q (z, P i). 



Thus for the elementary number density we simply have 



f;£Wl)=<Wx)(l + 0.1sm^), 



- A h ll 



(28) 



where d\l q {z,p\) is given by Eq. (21 1. For the toy 



model, we perform the MC simulations with only light 



z=0.6 




FIG. 7. The histograms (blue dots) for polarized number 
density D^o r v t (z, tp) for two values of z as a function of the 
azimuthal angle tp from NJL-jet framework using a toy model 
with a fixed number of hadrons emitted in each decay chain 
NLmks = 6. The minimum-x 2 fits with a functional form 




(a) 




D 



h/qt 



co + ci sin((y3) (red lines) are also depicted. 



(b) 



quarks in the quark-jet and pions in the final state, for 
simplicity. We perform several high-statistics simula- 
tions, where the number of produced hadrons in each 
decay chain, Nn n ). s , is fixed to a particular value. We 



use the relation ( 18 ) to extract the polarized number 
density D h / q ^(z, P±,(p) from the numbers of produced 
hadrons in intervals of the variables z, P\ and (p. For 
simplicity, in this section we will only consider the re- 
sults for the P^-integrated polarized quark fragmentation 
Dyr+fu^ ( 2 j </?) that can be obtained from Df l / q t(z, Pj_, (p) 
using the relation in Eq. Next, to extract the un- 
polarized and the 1/2 moment of the Collins functions 
Di(z) and iJT^ 1 (z), we perform a minimum-^ 2 fit to 
Df l / q j t(z,tp) using a form F(cg,ci) = cq + c\ sin(</?) for 
fixed values of z. This fitting method allows us to bet- 
ter account for the statistical fluctuations in the MC re- 
sults when extracting the Collins function. Notably, the 
fits describe the produced functions very well, yielding 
X over the number of degrees of freedom always in the 
vicinity of 1. A sample of such extraction for u — > tt° 
fragmentation is shown in Fig. [7J for simulations with 
Nhinks — 6. The histograms (blue dots) show the polar- 
ized number density / u t (z, </?) for two values of z as a 
function of the azimuthal angle (p. The minimum-^ 2 fits 
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(c) 



FIG. 8. Fitted values for unpolarized fragmentation function, 
D\, (a) and twice the Collins function 1/2 moment, 2H^ i - 1 ^ 2 \ 
(b) and their ratio (c) for and 7r versus z, produced 

by u quark from quark-jet framework using a toy model. The 
error bars show the uncertainties from the statistics and fits, 
and are only visible for the unfavored ratio at large values of 



with a functional form D h / q t — cq + c\ sin(y) (red lines) 
are also depicted. 

Then, using Dt(z) = 2nc and -2S q H^ {1/2) (z) = 
2itc\ [see Eq. ([Hjl], we plot the resulting fragmentation 
functions for u quark and Nu n k s — 6 as function of z 
in Fig. [HI setting, for example, S q = — 1. The plots in 
Fig. [8j a) depict the extracted, integrated fragmentation 
functions, which are in agreement with the results of the 
direct MC calculations in our previous work [TBTfTS] . 

The results for the moment of the Collins function, 
2i? 1 - L(1/2) , depicted in Fig. 8[b), exhibit very interesting 
features: both functions for u — > ir + and u — > 7r°, though 
peak at values of z ~ 0.6, decrease and u — > ir° changes 
the sign for small values of z. The unfavored function 
u — > 7r~ , generated solely by multihadron emission, has 
opposite sign and peak value at small z of comparable 
size to that of the favored ones at large z. A similar 
picture is found for the ratio 2H^' 2 ' /Di, depicted in 
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FIG. 9. Fitted values for 2# 1 L(1/2) (a) and the ratio 
2_ff 1 L ' 1/ ' 2 ' 1 / Di (b) for 7r + and tt~ versus z, produced by u 
quark from quark-jet framework using a toy model for in- 
creasing values of N L = N Links . 



Fig. ^Tc) , where notably the results for n° are below those 
for 7T^, while these ratios for the corresponding elemen- 
tary functions coincide. 

To better understand the results for the u — > tt + and 
u — > 7r~ Collins functions, in Fig. [9] we depict the results 
for N^in^s equal to 1, 2 and 6. The plots for 2if^ 1/ ' 2 ' in 
Fig. |9^a) show that the functions change very little for 
NLinks > 2 in the high z region. That is, the bulk of the 
model effects for z > 0.2 can be described with just two 
hadron emissions. We verified numerically with lower 
statistics runs, that the further increase in Nn n k s > 6 
affects the functions only at extremely low values of z, be- 
low the numerical discretization size of Az = 0.01 used 
in this study. A naive interpretation for the results is 
that the remnant quark has larger probability to have 
its spin antiparallel to that of the splitting quark, see 
Eq. (16), and typically has a small fraction of the initial 



light-cone momentum. Then in the emission step of the 
second hadron that affects the low z region and gener- 
ates the unfavored function, the angle if in the the Collins 
term in Eq. (27) in most cases will acquire an additional 



7r phase, yielding the results in Fig. [9] This picture re- 
peats for further hadron emissions, creating "destructive 
interference" in the small- z region, decreasing the Collins 
functions. The ratio 2 J ff 1 L(1/2) /-Di, depicted in Fig.^b), 
also exhibits the contrast in the low-z behavior of the 
unpolarized and Collins functions: while D±(z) grows 
roughly as 1/z, the Collins function oscillates to 0. It is 
worth mentioning that, for Nn^s = 1, the fitted value 
of the ratio in u — > ir + equals to the one set in the toy 
model as input to the MC. 



V. COLLINS FUNCTION FOR PIONS AND 
KAONS 

In this section we present the results of the full MC 
simulations with light and strange fragmenting quarks as 
well as pions and kaons as the produced hadrons. The ele- 
mentary splitting functions were taken from Eq. (27) and 
the simulations were done for values of Nn^s equal to 1, 
2 and 6. Again, it was checked with lower statistics runs, 
that the solutions are indistinguishable for N^^s > 6 
with the number of discretization points for z, P\ and tp 
used here. 



A. Pj integrated results 

We first present the results for the Pj-integrated 
Collins function, skipping the results for the unpolarized 
functions, D%(z), as they have been studied in detail in 
our previous work. 

The results for 2H 1 L(1/2) and the ratio 2iJ 1 L(1/2) /Z?i for 



the had rons produced by a u quark are shown in Figs. 10 
and [IT] respectively. The fitted values for 2H^ 1 ^ 2 ^ for 
pions show similar features to those from the toy model: 
the favored functions, positive and peaking at z ~ 0.65, 
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FIG. 10. Fitted values for 2H^ ( ~ 1 ^ 2 ^ for pions (a) and kaons 
(b), produced by u quark for Nunks = 6. 



FIG. 11. The ratio of the fitted values, 2H^ 1/2) /D 1 , for pions 
(a) and kaons (b), produced by u quark for N Links = 6. 



decrease and some change sign at lower values of z, os- 
cillating around zero as z — > 0. The unfavored function, 
u — > 7r~, is mostly negative and peaking at z ~ 0.2, with 
a peak value about one-third of that for u — >• tt + . The re- 
sults for kaons are similar to the those for pions, with the 
exception of u — > K + . This is a favored fragmentation, 
but the peak value of the unfavored function u — > K~ is 
only slightly less of that for the favored fragmentation. 
Also, since in our model charge and isospin symmetries 
are exact, the results for u — » K° coincide with the ones 
for u — > K~ and are omitted. 

Similar results for the s quark are shown in Fig. |12| 
The fitted values for 2iJ 1 L(1/2) for pions and unfavored 
kaon channels peak at value of z ~ 0.1 with values much 
smaller than those for the favored channels of K~ and 
K , which have a very broad peak and do not become 
negative at any value of z. 
The dependence on 



Links 



can be seen in Fig. 13 



Again, as for the toy model, we notice that the results 
change very little for z > 0.2 for Nn n k s ^ 2. The 
increase in Ni,inka affects the functions at decreasingly 



smaller values of z. 

B. TMD of Collins function 

In this section we present the full results for the TMD 
Collins function. As the z dependence has been stud- 
ied in detail in the previous sections, here we will only 
show the P\ dependence of the Collins function for sev- 
eral fixed values of z. We noticed that the polar angle 
dependencies of the polarized number densities are not 
very accurately described by just the sine modulation 
term corresponding to the Collins term for some values 
of z and Fj\ I n our calculations with Nn n k s = 6, we 
found that a fourth order polynomial in sin(y) 



D h/qt {z, Pi,ip) = J2 c n(z, Pi) sin" p (29) 

is sufficient to achieve fits for all slices with fixed z, P\ 
with x 2 per degree of freedom in the vicinity of 1. Here we 
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FIG. 12. Fitted values for 2# 1 L(1/2) (a) and 2# 1 L(1/2) /-Di (b), 
for hadrons produced by s quark for N^i n k a — 6. 



FIG. 14. The unpolarized fragmentation function vs P ± for 
pions produced by a u quark with z = 0.7 (a) and z = 0.1 (b) 
for N Linka = 6. 




FIG. 13. Fitted values for 2# 1 ±(1/2) for K+ and K~ versus 
z, produced by u for increasing values of Nl = N Links . 



identify cq with the unpolarized and c\ with the Collins 
term of Eq. The higher powers of the sin(ip) in the 
polarized quark fragmentation functions are induced by 
the multiple hadron emissions. A toy model analysis 
showed that there are two sources for these effects: the 
remnant quark transverse momentum distribution modu- 
lation from the recoil in the elementary hadron emissions, 
and an additional modulation of the same quark distri- 
bution due to the ip dependence of the quark spin flip 
probabilities. The detailed description and explanation 
of these additional terms are presented in our forthcom- 
ing publications [221 130], where we show that these effects 
are a genuine feature of the quark-jet model and do not 
depend on the particular form of the elementary polar- 
ized fragmentation function. The plots in Figs. [T4j [15] 
and [16] depict the results for unpolarized, Collins frag- 
mentation functions and their ratios respectively for pi- 
ons produced by an initial u quark with two fixed values 
of z equal to 0.7 and 0.1. For z = 0.7, the results are 
mostly affected by a single hadron emission, thus the P]_ 
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FIG. 15. The Collins fragmentation function vs P]_ for pions 
produced by a u quark with z = 0.7 (a) and z — 0.1 (b) for 

Nhinks = 6. 



FIG. 16. The ratio of Collins fragmentation function to un- 
polarized fragmentation function vs P]_ for pions produced by 
a u quark with z = 0.7 (a) and z = 0.1 (b) for NLinks = 6. 



dependence of the Collins function is somewhat similar 
to one of the unpolarized fragmentations (which, in turn 
are reasonably well described via Gaussian function for 
small values of at a fixed z, see p]), but peak at bein S smaller - 
a small nonzero value of Pj. The results are quite dif- 
ferent for z — 0.1, where the unpolarized functions for 
all the pions are roughly equal, while the Collins func- 
tion for all the pions have similar shapes but differ in 
magnitude, decreasing to a negative value and peaking 
at Pj_ ~ O.lGeV 2 , then start to increase and become 
positive for small values of P]_. Thus, it follows from 
our model that while the transverse momentum depen- 
dence of the unpolarized fragmentation functions can be 
to a good approximation described by a Gaussian func- 
tion with z-dependent width, the Collins functions have 
much different shapes. 



The plots in Figs. [T7j [18| and [19] depict the analogous 
results, but for kaons produced by an initial u quark. 
The results are similar to those for the pions, except for 
the magnitudes of the unpolarized and Collins functions 



The plots in Figs. |20| [21] and [22] depict the results for 
unpolarized, Collins fragmentation functions and their 
ratios respectively for pions and kaons produced by an 
initial s quark with two fixed values of z equal to 0.7 and 
0.1. In our model the isospin symmetry is considered to 
be exact, thus the calculated values for all the pions are 
equal to each other, thus only 7r + is depicted. Similarly, 
the results for the kaons that can be transformed to each 
other by an isospin rotation coincide and are omitted as 
well. The results are similar to those for the u quark, 
except for the strong dominance of the favored fragmen- 
tation functions at z — 0.7. At z = 0.1, where multiple 
hadron emission contributions dominate, the unfavored 
fragmentation functions become larger than the favored 
ones, in particular for the production of pions. 



12 




0.4 0.6 

P[ (GeV 2 ) 

(a) 



0.2 



P[ (GeV 2 ) 

(a) 



u->h 

K°,z=0.7 

K+,z=0.7 

K\ z=0.7 



0.4 0.6 

Pi (GeV 2 ) 

(b) 



K",z=0.1 
K+,z=0.1 
K,z=O.I 




().]() 



-0.05 



1 1 1 1 1 1 1 

] 

■ 


1 1 1 
u->h 


K°,z=0.1 - 


1 
\ 


K+,z=0.1 - 


ST, z=0.l ' 












X 






1.1,1 


l 



0.4 0.6 

P[ (GeV 2 ) 

00 



FIG. 17. The unpolarized fragmentation function vs P ± for 
kaons produced by a u quark with z = 0.7 (a) and z = 0.1 
(b) for N Link s = 6. 



FIG. 18. The Collins fragmentation function vs P ± for kaons 
produced by a u quark with z = 0.7 (a) and z = 0.1 (b) for 

NLinks = 6. 



C. The Schafer-Teryaev sum rule 

The Schafer-Teryaev sum rule was originally proposed 
in Ref. [51] . where the authors used the simple argu- 
ments of the transverse momentum conservation in the 
hadronization process to derive 



h 



(30) 



where the first moment of the Collins function is defined 
as 



H 



■l(i) 

l,(h/q) 



dPl 



2zrrih 



H 



±h/q 



(31) 



Later, in Ref. [32], it was proven explicitly using 
the quark correlator functions, along with the quark 
to hadron Collins function moments that the quark-to- 
quark Collins functions should also be included for the 



sum rule to be satisfied 



dz2H^ /q) (z) = 0. (32) 



Our model results both for the toy model and the full 
calculation, show that the naive Schafer-Teryaev sum 
rule cannot be satisfied, even though the transverse mo- 
mentum conservation is explicitly satisfied in our simu- 
lations. For example, the sum of the terms in Eq. (30) 



for u quark is ST U — 0.07, while for s quark ST S = 0.21. 

A simple explanation for this is that with multiple 
hadron emissions the quark quickly loses almost all of its 
initial light-cone momentum, but gains a nonzero aver- 
age transverse momentum from the recoil of the emitted 
hadrons. This can be easily seen from the fact that the 
solutions of unpolarized integrated fragmentation func- 
tions, -Di(z), change only in the extremely small region 
of z, after only a few hadron emissions, as described in 
Ref. [T7] . On the other hand, the average transverse mo- 
mentum of the emitted hadrons tends to level off at a 
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FIG. 19. The ratio of Collins fragmentation function to un- 
polarized fragmentation function vs Pj_ for kaons produced 
byan quark with z = 0.7 (a) and z = 0.1 (b) for Nhinks = 6. 

certain nonzero value as z — > 0, as shown in Fig. 14 of 
[18j , hinting at a similar value for the average transverse 
momentum of the remnant quark. In the NJL-jet model 
we assume that the remnant quark hadronizes with the 
slow, colored fragments of the initial struck hadron in 
SIDIS, or the remnant of the antiquark in e + e~ reac- 
tions, with the product unobserved because of the small 
values of remaining z it carries. Thus, while the effects of 
the remnant quark after multiple hadron emissions can 
be neglected when considering the longitudinal momen- 
tum sum rules for the unpolarized fragmentations, they 
are essential for the transverse momentum sum rule for 
the Collins function, as shown in |32j . 

VI. CONCLUSIONS 

In this article we calculated the Collins fragmentation 
function in the NJL-jet model. This was accomplished 
by extending the model to include the transverse polar- 
ization of the fragmenting quarks and calculating the 
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FIG. 20. The unpolarized fragmentation function vs Pi for 
pions and kaons produced by a s quark with z — 0.7 (a) and 
z = 0.1 (b) for N Links = 6. 



probability of the quark spin flip at each hadron emis- 
sion. Then, we extended our Monte Carlo framework 
to accommodate for the spin of the quarks in the jet, 
and polar angle of the produced hadrons' transverse mo- 
menta with respect to the direction of the initial quark's 
spin. The polar-angle dependent elementary polarized 
fragmentations functions, that are a sum of elementary 
unpolarized and Collins functions, were used as inputs for 
the MC simulations. Here the unpolarized function was 
taken the one used earlier in NJL-jet model, while the 
Collins function was taken from calculations within the 
spectator model of |11) , with the parameters and regular- 
ization of transverse momentum integrals of the NJL-jet 
model. Using this input, we calculated the correspond- 
ing number densities for the produced hadrons in MC 
simulations, where we fixed the number of the produced 
hadrons in each quark-jet (N L i n k s )- Then, we used the 
form of the polar angle dependence of the full number 
densities of Eq. ^ to separate the unpolarized and the 
Collins functions. 
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FIG. 21. The Collins fragmentation function vs P ± for pions 
and kaons produced by a s quark with z = 0.7 (a) and z = 0.1 
(b) for N Link s = 6. 



FIG. 22. The ratio of Collins fragmentation function to un- 
polarized fragmentation vs P\ for pions and kaons produced 
by a s quark with z — 0.7 (a) and z — 0.1 (b) for NLinks ~ 6. 



Results for two scenarios were presented. First, we 
showed the results for Pj_ integrated case using a toy 
model, where we assumed the form of the elementary 
Collins term to be that of the unpolarized one times 
0.1. Also, the simulations were performed with only light 
quarks and pions for simplicity. The results for the so- 
lutions of Hi^~' 2 ' (z) showed that a nonzero unfavored 
Collins function can be generated in our model, with 
a sign opposite to that of the favored one. Also, both 
the unfavored and favored Collins functions oscillate for 
the lower values of z. Second, we presented the results 
with the full model calculations of both unpolarized and 
Collins functions with light and strange quarks, as well as 
pions and kaons. The resulting full solutions of the 1/2 
moment of the Collins functions exhibit features simi- 
lar to the ones from the toy model. Interestingly, for 
the solutions for the kaons produced by a u quark, the 
magnitudes of the unfavored solutions were only slightly 
smaller than the favored ones. 

Finally, we studied the P\ dependence of the Collins 



function for several fixed values of z. Here we noticed 
that the multiple hadron emissions produce modulations 
of the polarized number densities, which contain several 
powers of sin(y). These higher order modulations are 
a genuine feature of the quark-jet hadronization mecha- 
nism and neither depend on the particular forms for the 
elementary fragmentation functions nor are they artifacts 
of our Monte Carlo simulations. A detailed study, which 
will be presented in our forthcoming publications [251150] . 
show that this effect can be completely described by the 
azimuthal modulation of the transverse momentum dis- 
tribution of the remnant quark in each elementary frag- 
mentation process and to the azimuthal modulations aris- 
ing from the quark spin flip probability. The fitting func- 
tions for obtaining the unpolarized and Collins functions, 



described in Sec. IV were modified using Eq. ( 29 ) to ac 



re unmte- 



count for this effect, which shows up only in t 
grated TMD fragmentation functions. Our results show 
that the Gaussian function cannot be used to reliably 
model the Collins functions in any regions of z. 
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We also investigated the Schafer-Teryaev sum rule for 
our results, and found that the naive sum rule of Ref. |31j 
cannot be satisfied with our solutions, even though the 
transverse momentum conservation is explicitly enforced 
in our MC framework. The omission of the remnant 
quark transverse momentum is the key here, as it was 
proven explicitly from the definition of the Collins func- 
tion in Ref. |32j that the sum rule holds if the quark- 
to-quark Collins functions are also accounted for. We 
argued that the remnant quark in the hadronization pro- 
cess after several hadron emissions, though carrying a 
minuscule fraction of the initial quark's light-cone mo- 
mentum, acquires a significant transverse momentum on 
average, thus needs to be included in this sum rule that 
manifests the transverse momentum conservation in the 
hadronization process. 

Our results for Collins functions have some distinctive 
features that were approximately observed in the exper- 
iment: the similar size and the opposite sign for the 1/2 



moments of favored and unfavored ones. At this stage, 
however, a direct comparison with the experimental data 
is not possible because we did not perform the QCD evo- 
lution on our results, as these equations for the Collins 
function are not yet known. The recent work on the 
evolution of TMD distribution and fragmentation func- 
tions of Ref. |33| pave the way for this, and we plan to 
implement the QCD evolution of Collins functions in our 
forthcoming work as the relevant formalism is developed. 
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